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(57) ABSTRACT 

A polarization imaging apparatus measures the Stokes image 
of a sample. The apparatus consists of an optical lens set, a 
first variable phase retarder (VPR) with its optical axis 
aligned 22. 5°, a second variable phase retarder with its optical 
axis aligned 45°, a linear polarizer, a imaging sensor for 
sensing the intensity images of the sample, a controller and a 
computer. Two variable phase retarders were controlled inde- 
pendently by a computer through a controller unit which 
generates a sequential of voltages to control the phase retar- 
dations of the first and second variable phase retarders. A 
auto-calibration procedure was incorporated into the polar- 
ization imaging apparatus to correct the misalignment of first 
and second VPRs, as well as the half-wave voltage of the 
VPRs. A set of four intensity images, I 0 , 1 15 1 2 and I 3 of the 
sample were captured by imaging sensor when the phase 
retardations ofVPRs were set at (0,0), (it,0), (it, it) and (it/2,it), 
respectively. Then four Stokes components of a Stokes image, 
S 0 , Sj, S 2 and S 3 were calculated using the four intensity 
images. 

12 Claims, 5 Drawing Sheets 
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POLARIZATION IMAGING APPARATUS 
WITH AUTO-C ALIBRATION 

CROSS REFERENCE TO RELATED 
APPLICATION 

This is a continuation-in-part of U.S. patent application 
Ser. No. 1 1/701,738, filed on Feb. 2, 2007 now U.S. Pat. No. 
7,701,561, which claims the benefit of provisional applica- 
tion Ser. No. 60/772,779 filed on Feb. 13, 2006, the entire 
teachings of all of which are incorporated herein by this 
reference. 

STATEMENT OF GOVERNMENTAL INTEREST 

This invention was made partially with Government sup- 
ports under grant No. NNJ05JC13C awarded by NASA and 
No. DE-FG02-08ER85109 awarded by Department of 
Energy. The government has certain rights in this invention. 

BACKGROUND OF THE INVENTION 

1. Technical Field 

The invention relates to optical imaging apparatus, particu- 
larly related to polarization imaging apparatus, more specifi- 
cally related to a high speed polarization imaging system that 
can measure Stokes components from a sample at different 
optical wavelengths. 

2. Technical Background 

The light scattered by a tissue has interacted with the ultra- 
structure of the tissue, which imprinted some intrinsic prop- 
erties of the tissue. Tissue ultrastructure extends from mem- 
branes to membrane aggregates to collagen fibers to nuclei to 
cells. Photons are most strongly scattered by those structures 
whose size matches the photon wavelength. It has been dem- 
onstrated that light scattering can provide structural and func- 
tional information about the tissue. One important biomedi- 
cal application of optical imaging and spectroscopy is non- 
invasive or minimally invasive detection of pre-cancerous and 
early cancerous changes in human epithelium, such as dys- 
plasia or carcinoma in situ. 

In recent years there has been an increasing interest in the 
propagation of polarized light in randomly scattering media. 
The investigation of backscattered light is of particular inter- 
est since most medical applications aimed at the in-vivo char- 
acterization of biological tissue rely on backscattered light. 
By recording the spatially dependent response of a medium to 
a polarized point source, one may obtain infomiation about 
the scattering particles that are not accessible to mere inten- 
sity measurements. A diagnostic imaging modality based on 
near-infrared (NIR) radiation offers several potential advan- 
tages over existing radiological techniques. First, the radia- 
tion is non-ionizing, and therefore reasonable doses can be 
repeatedly employed without harm to the patient. Second, 
optical methods offer the potential to differentiate between 
soft tissues, due to their different absorption or scatter at NIR 
wavelengths that are indistinguishable using other modali- 
ties. And third, specific absorption by natural chromophores 
(such as oxy-haemoglobin) allows functional information to 
be obtained. NIR imaging research has focused on a variety of 
possible clinical applications. 

The fact that the polarization state of the light contains 
useful information has been shown in many literatures, for 
example. Rahmann and Canterakis in “Reconstruction of 
specular surfaces using polarization imaging,” IEEE Com- 
puter Society Conference on Computer Vision and Pattern 
Recognition , vol. 1, pp. 149-155, 2001, describe how the 
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polarization state of light canbe used forthe reconstruction of 
specular surfaces to determine the shape of 3-D objects. They 
use the fact that light reflected by dielectrics and metals 
becomes partially linearly polarized and that the direction of 
5 polarization depends on the orientation of the reflecting sur- 
face. Demos and Alfano in “Deep subsurface imaging in 
tissues using spectral and polarization filtering,” Optics 
Express, vol. 7, no. 1, pp. 23-28, 2000, demonstrate a tech- 
nique based on polarization imaging that allows for optical 
to imaging of the surface as well as structures beneath the sur- 
face. 

There have been several polarization imaging schemes 
published. In U.S. Pat. No. 6,437,856, inventor Steven Louis 
Jacques, which issued Aug. 20, 2002, and U.S. Pat. No. 5,847, 
15 394, inventors Robert R. Alfano, et al., which issued Dec. 8, 
1 998, a set of measurements at different polarization orien- 
tations are taken to render a new image that is independent of 
the light reflected from the surface of a tissue sample and that 
is dependent of the light scattered from deep tissue layers. 
20 Especially, a linearly parallel polarized light is used for illu- 
mination and two images are acquired, one image selecting 
linearly parallel (Par) polarized light (i.e., parallel to the light 
source-tissue-camera plane) and one image selecting linearly 
perpendicular (Per) polarized light (i.e., perpendicular to the 
25 fight source-tissue-camera plane). A new image is obtained 
by (Par-Per) or (Par-Per)/(Par+Per). 

Many polarimetric sensing technologies have been devel- 
oped to capture the Stokes polarization information. People 
use rotating retarder, rotating polarizer to obtain Stokes 
30 parameters from several successive frames of image. In these 
mechanically-rotation approaches, however, the exact spatial 
registration between various frames is difficult due to changes 
of viewing angle and the polarization properties during the 
period of between successive frames in most target detection 
35 applications. 

In U.S. Pat. No. 6,798,514, inventor James Maurice 
Daniels, which issued Sep. 28, 2004, the light passes through 
two liquid crystal waveplates and a polarizing filter before 
falling on a light sensitive element device to measure intensity 
40 of light in four different polarizations. However, the same 
difficulty also applied to the approach using liquid crystal for 
electro-optic polarization modulation, which usually has 
response time in the order of 100 ms. Another drawback for 
the liquid crystal approach is that polarization-rotating liquid 
45 crystals are not commercially available in many important IR 
bands. 

SUMMARY OF THE INVENTION 

50 It is therefore the main objects of the invention to provide 
a polarization imaging apparatus which can consistently pro- 
vide high speed and broadband Stokes images in real-time. 

It is also the object of the invention to provide a polariza- 
tion imaging apparatus with an auto-calibration procedure. 
55 The key feature of the polarization imaging apparatus is 
that it provides fast adjustable optical phase control through 
high speed electro-optical variable phase retarders made from 
transparent electro-optic ceramics. 

A preferred embodiment of the polarization imaging appa- 
60 ratus is shown in FIG. 1, which consists of a polarization 
imaging sensing unit 100, an electronic driver 101 and a 
computer 102 . The imaging sensing unit 100 consists of an 
optical assembly 1 1 , a linear polarizer 14 with its optical axis 
18 , a first variable phase retarder (VPR) 12 with its optical 
65 axis 16 aligned 22.5° to axis 18 , a second variable phase 
retarder 13 with its optical axis 17 aligned 45° to axis 18 , and 
a imaging sensor 15 for sensing the intensity images of the 
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sample. The imaging sensor 15 can be a focal plan array 
(FPA) detector, a CMOS or a CCD sensor. Two VPRs 12 and 
13 were controlled independently by a computer 102 through 
a controller unit 101 which generates a sequential of voltages 
to control the phase retardations of VPRs 12 and 13 . A set of 
four intensity images, I 0 , I x , I 2 and I 3 of the sample were 
captured by imaging sensor 15 when the phase retardations of 
VPRs 12 and 13 were set at (0,0), (jt,0), and (it/2,it), 
respectively Then four Stokes components of a Stokes image, 
S 0 , S l5 S 2 and S 3 can be calculated by computer 102 using the 
following formula: 

5 0 = 0.5 x (/<, + /,) (1) 

51 =0.5 x(/ 0 -/i) 

' S 2 = 0.5x(I 2 -S o ) 

. s 3 = 0.5 x(/ 3 - So) 

Another key feature of the invention is that a simplified 
time sequential control and measurement algorithm which 
provides a fast measurement and converts sequential intensity 
images to Stokes images. 

In one aspect, the present invention has a calibration pro- 
cedure that needs one known polarization state to accurately 
calculate all of the parameters, such as the orientation of first 
and second variable phase retarter al and a2, respectively, as 
well as half-wave phase retardance 61” of first variable phase 
retarder, quarter- wave phase retardance 6 1’"' 2 of first variable 
phase retarder, and half-wave phase retardance 52™ of second 
variable phase retarder. 

In one aspect, the present invention includes the ferro- 
electric materials, such as PLZT, La modified PMN-PT, and 
La modified PZN-PT to form the electric voltage controlled 
ferro-electric variable phase retarder. 

Referring to FIG. 5 is another preferred embodiment of the 
polarization imaging apparatus consists of an illumination 
unit 50 , a polarization imaging sensing unit 500 , a controller 
unit 101 and a computer 102 . The said illumination unit and 
imaging sensing unit are aligned in a reflection configuration. 
In addition to a light source, the illuminator may consist of 
parts such as a collimator, a polarizer or polarization modu- 
lator, and a fixed wavelength filter or tunable wavelength 
filter. Tunable optical band-pass filter is used to produce 
desired wavelength for the imaging, from visible to infrared 
(IR). The illuminator can produce a linear or arbitrary polar- 
ized light beam. 

The illumination unit 50 , and a polarization imaging sens- 
ing unit 500 in FIG. 5 can be further aligned in a transmission 
configuration, as illustrated in FIG. 6. 

Additional features and advantages of the invention will be 
set forth in the detailed description which follows, and in part 
will be readily apparent to those skilled in the art from that 
description or recognized by practicing the invention as 
described herein, including the detailed description which 
follows, the claims, as well as the appended drawings. 

It is to be understood that both the foregoing general 
description and the following detailed description present 
embodiments of the invention, and are intended to provide an 
overview or framework for understanding the nature and 
character of the invention as it is claimed. The accompanying 
drawings are included to provide a further understanding of 
the invention, and are incorporated into and constitute a part 
of this specification. The drawings illustrate various embodi- 
ments of the invention and together with the description serve 
to explain the principles and operations of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is the schematic of a preferred embodiment of 
polarization imaging apparatus with variable retarders in 
accordance with the present invention; 

FIG. 2 shows the measurement steps for polarization imag- 
ing apparatus in accordance with the present invention; 

FIG. 3 shows the schematic drawing of a high speed vari- 
111 able retarder using a electro-optic ceramic; 

FIG. 4 shows measured responses of a high speed variable 
retarder; 

FIG. 5 is a schematic view of a preferred embodiment of 
15 reflective polarization imaging apparatus in accordance with 
the present invention; 

FIG. 6 is a schematic view of a preferred embodiment of 
transmissive polarization imaging apparatus in accordance 
with the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

It is to be understood that the invention is not limited in its 
application to the details of the particular arrangement shown 
since the invention is capable of other embodiments. Also, the 
terminology used herein is for the purpose of description and 
30 not of limitation. 

FIG. 1 shows a preferred embodiment of the polarization 
imaging apparatus in accordance with the present invention. 
It consists of an imaging sensing unit 1 00 , an electronic driver 

101 and a computer 102 . The imaging sensing unit 100 fur- 
35 ther consists of an optical assembly 11 , a polarizer 14 with its 

optical axis 18 , a first variable phase retarder (VPR) 12 with 
its optical axis 16 aligned 22.5° to axis 18 , a second variable 
phase retarder 13 with its optical axis 17 aligned 45° to axis 
40 18 , a imaging sensor 15 for sensing the intensity images of the 
sample. The imaging sensor 15 can be a focal plan array 
(FPA) detector, a CMOS or a CCD sensor. The imaging 
sensing unit then is assembled inside a metal host (not 
shown). Two VPRs 12 and 13 were controlled independently 
45 by a computer 102 through a controller 101 which generates 
a sequential of voltages to control the phase retardations of 
the VPRs 12 and 13 . The intensity images were captured 
through the imaging sensor 15 , and processed by computer 

102 to convert the sequential intensity images to Stokes 
50 images. 

Using Mueller’s calculus, we can obtain the Stokes param- 
eters of an input light 10 by measuring the output light inten- 
sity of the sensor 15 . The transformation done by the phase 
retarders can be described in the following equation, 

^r(So,^i,„»&, m &,„J r =K 2 *M rl ( 2 ) 

where S,„=(S 0 Si S 2 S 3 ) r is the unknown input Stokes vector 
of the target 10 , which consists of four Stokes components. 
60 The parameter S 0 is the intensity of the light (including un- 
polarized component), Sj is the intensity difference between 
horizontal linearly polarized and the vertical linearly polar- 
ized components, S 2 is the intensity difference between 45° 
65 linearly polarized and the 135° linearly polarized compo- 
nents, and S 3 is the difference between left- and right-circu- 
larly polarized intensities. 
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The quantities M rl , and M,, 2 are the Mueller matrices rep- 
resenting two variable phase retarders 12 and 13 , respectively. 
These are given for the specific orientation by 


1 

0 


M rl m = 


0 


1 


M^si) = - 


0 0 
i(l+cos(<51)) i(l-cos(tfl)) 

i(l -cos(<51)) i(l+cos(tfl)) 

1 1 
— — sin(<51) — — sin(<51) 

V2 V2 

1 0 0 0 ' 

0 cos((52) 0 — sin(<52) 

0 0 10 
0 sin(<52) 0 cos(<52) 


0 

1 

— sin((51) 

V2 

1 

— — sin(<51) 


cos(<51) 


5 

(3) 


10 


(4) 

15 


where 61 and 62 are the two phase retardance introduced 20 
by the two variable phase retarders 12 and 13, respectively. 

To be able to measure the four components of the input 
Stokes vector, the variable phase retarders need to be set to at 
least four different states labeled as k. In the following, the 
linkage between input state Stokes parameters and measured 23 
light intensity is detailed. The Mueller matrix for the cas- 
caded variable retarders can be generally cast into an alterna- 
tive form so that the state of polarization incident on the 
polarizer can be expressed as, 


' S(),out 


S(), in 


’moo 

0 

0 

0 ' 

' So,in ' 

Si ,out 

= M ■ 

Si, in 
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S 2 ,out 
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m 2 i 

m 22 

m 2 3 

S 2 ,in 

, ^3, out , 


. S 3 , in . 


0 

m 3l 

m 32 

m 33 , 

,S3,in , 


The actually detected signal through the polarizer is essen- 
tially 


45 

Without worrying about a global factor, the factor Vi was 
thrown away. Thus the intensity can be expressed in terms of 
three of the Muller matrix elements, 

,S Ui „+m 12 S 2 i „+m I 3 S 3J „ (6) 50 

It is clear that the output intensity can be linked to the input 
Stokes parameters through generating four known states of 
polarization at the input so that we have, 


'h' 
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where mT is the Mueller matrix element at the k-th state of 
variable retarders, and I, is the intensity measured with the i-th 
state of variable retarder. From Equation (7), it can be noticed 65 
that once the matrix A is known, the input Stokes parameters 
can be found out the inverse the matrix so that, 



(8) 


The very question now is to find out the matrix A. It can be 
calculated from each Muller matrix in the system or, it can be 
done by generating four known states of polarization so that, 


r>l r >2 c <3 rr 4 

°0,m °0,/n °0 ,in °0 ,in 


I 1 i 2 / 3 I 4 

l l l l l l *1 

Sljn S?,„ sfj„ sf j„ 


/ 2 / 4 

l 2 l 2 l 2 ‘2 

ril ri 2 r <3 rr 4 

°2 ,in °2,in °2 ,in °2 ,in 


/I / 2 / 3 I 4 

/ 3 l 3 l 3 h 

ol o 2 rt 3 o 4 

03 j n Oj j n A3 [ n 03 


,1 ,2 ,3 ,4 

‘4 l 4 l 4 ‘4 


where S, J is the i-th Stokes parameter generate by the j-th 
retardance. Thus we are able to link the input state polariza- 
tion to the intensity measured at four retardance states. The 
Mueller matrix elements needed for finding the input state of 
polarization can be figured out through the use of four known 
state of polarization. 

The measurement procedure is based on specifying above 
general principle for practical applications. As the configura- 
tion shown in FIG. 1, we can find out the useful Muller matrix 
elements (according to equation (6)) as follows: 

»'oo =1 


m, j=cos 2 (6l/2)*cos(62)-sin(61)-sin(62)72 


m i 2 =sin 2 ( 6 l 2 ) 


m i3 =-cos 2 (6l/2)*sin(62)-sin(6l)-cos(82)/V f 2 (10) 

where 51 and 62 are the retardance of two VPRs 12 and 13, 
respectively. Although their values could be many different 
combinations for finding out the Stokes vector, the following 
4-state procedure is considered to be the simplest yet most 
effective approach. 

Referring to FIG. 2 shows the procedure diagram to 
accomplish one measurement of Stokes parameters. 

Four intensity images (I 0 -I 3 ) are needed and taken sequen- 
tially, with two variable phase retarders 12 and 13 controlled 
by a sequence of voltages which were generated by controller 
101 . Phase retardation changes 21 and 22 ofVPR 12 and VPR 
13, respectively. I 0 is taken with 0 volt applied to both 12 and 
13. Ij is taken with the half-wave voltage V„ applied to 12 and 
0 volt applied to 13. I 2 is taken with V, is applied to both 
retarder 12 and 13. I 3 is taken with a quarter-wave voltage 
V* /2 applied to 12 and applied to 13. The four Stokes 
vector images S 0 -S 3 can be calculated from the intensity as 
follows: 


5b=0.5x(/o + /i) (11) 

5) = 0.5 x (/<,-/,) 

5 2 = 0.5 x (/ 2 - 80 ) 

5 3 = 0.5 x (/ 3 — 80 ) 


Table 1 listed the state of retardance of first and second 
variable phase retarders and the corresponding four intensi- 
ties measured in each equal length time slot. The intensities 
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also can be expressed as a function of Stokes parameters 
S 0 -S 3 of the input Stokes vector according to equation (2). 
Table 1 also listed the formula for extracting the Stokes 
parameters S 0 -S 3 from these intensities. 


TABLE 1 


Time slot 

1 

2 

3 

4 

State ofVPRl ( 61 ) 

0 

jt 

jt 

jt/2 

State of VPR 2 ( 62 ) 

0 

0 

jt 

jt 

Measured intensity 

Io 

I x 

h 


Relation to Stokes 

So./„ + S w „ So,*- 

Si j n ^0,in + ^2Jn 


Derived parameters 


'S 0 : 

= 0 . 5 x(I 0 + Ii) 




Si : 

= 0 . 5 x(I 0 -Ii) 




s 2 = 

= 0.5 x (I 2 - So) 




s 3 = 

= 0.5 x (I3 — Sq) 



known polarization incident light, such as a horizontal polar- 
ized light: 


5 

Sout{7 r/2,0) 



Sout(7T, 0) 



SoUt,(0,7[) 
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^OUt(77,77) 


10 

. SoUt(7t/ 2jt) . 



where 

15 

A,;„ 

S 2 .in 

A,i„ 

20 


, /i/2,0 /i/2,0 ii/2,0 

1 mil m 12 '«13 

| 7T,0 77,0 7T,0 

1 »iu m, 2 nii 3 

, O^r O^r 0,7r 

1 in n m 12 mi 3 

t 71,71 71,71 71,77 

1 nin mil m3 

1 m ^/ 2 ^ m 7r / 2 ’ 7r 

i m n m 12 m 13 


5o,i 

S'l./n 

ft* 


(15) 


With this simple converting formula, Stokes images can be 
obtained using very simple calculation from the intensity 
images captured by an imaging sensor. 

Considering the misalignment of both retarders VPR 12 
and VPR 13, and the inaccuracy of their actual half-wave 
voltages in a polarization system, an auto-calibration tech- 
nique is necessary to maintain its accuracy. The auto-calibra- 
tion procedures are for correcting the errors caused by mis- 
alignment the optical axis 1 6 ofVPR 12 and optical axis 17 of 
VPR 13, and the inaccuracy of the voltages applied to each of 
VPRs. 

Assuming: 

al the actual orientation of the optical axis 16 of variable 
phase retarder VPR 12. The ideal value is 22.5°; 

a2 the actual orientation of the optical axis 17 of variable 
phase retarder VPR 13. The ideal value is 45°; 

8 1 51 the actual phase retardance when a half-wave voltage 
V, is applied to VPR 12. 

The ideal phase retardance is Jt; 

8 1 31/2 the actual phase retardance when a quarter-wave volt- 
age V x/2 is applied to VPR 12. The ideal phase retar- 
dance is n/2; and 

82 31 the actual phase retardance when a half-wave voltage 
V, is applied to VPR 13 . The ideal phase retardance is it. 

where initial quarter-wave voltage and half-ware voltages 
ofVPR 12 and VPR 13 are measured before assembly. 

Hie components of Mueller matrix in Equation (5) will be: 

m,,=(cos 2 +2al+sm 2 2al cos 61 )(cos 2 2 a 2 +sin 2 2 a 2 
cos 62 )+sin 2al cos 2al(l-cos Sl)sin 2a2 cos 
2a2(l-cos 62 )-sin 2al sin 61 sin 2a2 sin 62 (12) 

Wi2 = ( c °s 2 2al+sin 2 2al cos 6l)sin 2a2 cos 2a2(l- 
cos 62 )+sin 2al cos 2al(l-cos 6l)(sin 2 2a2+ 
cos 2 2a2 cos 62 )+sin 2al sin 61 cos 2a2 sin 62 ( 13 ) 

w 13 =-(cos 2 2al+sin 2 2al cos 6l)sin 2a2 sin 62 +sin 
2al cos 2al (1-cos 6l)cos 2a2 sin 62 -sin 2al 
sin 61 cos 62 ( 14 ) 

The parameter al, a2, SI 3 ® 72 , 61 31 and 62 31 can be solved 
using a following sets of measurement procedure with a 


is the Stokes matrix of the input light and 


25 $out(] r/2,0) 

Sout(7C,0) 

Sout,(0jr) 

^out(n,7T) 

. S 0 ut(7 r/2^r) . 

30 

are measured light output. The Mueller matrix components 
are, according to equations (12) to (13): 

(1 ) m,, 31 2 ’°, m 12 ™ /2,0 , m^" 72 ’ 0 represent the m values while a 
quarter-wave voltage and a zero are applied to VPR 12 and 

35 VPR 13, respectively. The detected output is S out ( X / 2 ,oy 

(2) m,, 31 ’ 0 , m 12 ’ I, °, m 13 ™’° represent the m values while a 
half-wave voltage and a zero voltage are applied to VPR 12 
and VPR 13, respectively. The detected output is S OIrf(3t0 ). 

(3) m, ,°’ 3 ', 111 , 2 °’’', m^ 0 ’ 31 represent the m values while a zero 
40 and a half-wave voltage are applied to VPR 12 and VPR 13, 

respectively. The detected output is S o „, (0 „- ) 

(4) , m, 111 , 2 ’'’’', m,/ 1 ’ 31 represent the 111 values while half- 

wave voltages are applied to both VPR 12 and VPR 13. The 
detected output is S oul(x ^ 

45 (5) m, | , m^” , m, 3 2,3t represent the 111 values while a 

quarter-wave voltage and a half-wave voltage are applied to 
VPR 12 and VPR 13, respectively. The detected output is 
q 

out(jtjt/2y 

The Stokes vector of an unknown incident image (S 0 , S, , 
so S 2 , S 3 f can be obtained using the equation (16): 


So 
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where S ou2(0j0) , S out(nfi) , S out( „^ and S 0 „, (3t/2j3t) outputs 
60 when the first VPR 12 and second VPR 13 have the phase 
retardation of (0,0), (jt,0), (it, it), and (jt/2,jt), respectively; 
and 111 values can be calculated according to equations (1 2 ) to 
(13) using calibrated VPR parameters al, a2, SI 31 , 8U 2 and 
82V 

65 It should be pointed out that the switch of the relative 
positions ofVPR 12 and VPR 13 will not affect the outcome 
of the Stokes images according to equation (2). 
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It is to be appreciated that the orientation configuration 
between the two variable phase plates 12 and 13 might be 
achieved with an intervening fixed waveplate, so that physical 
orientation of the two plates could be less constrained. Hence, 
the 22.5 degree orientation of the second variable phase 5 
retarder 13 with respect to the first variable phase retarder 12 
is achieved by any combination of relative physical rotation 
of the second variable phase retarder 13 with respect to the 
first retarder 12 and rotation of the polarization reference 
frame between the two variable phase retarders 12 and 13 by 1 0 
a fixed phase retarding plate (not shown but included as part 
of the present invention). 

Any other orientation configuration of first variable phase 
retarder 12 and second variable phase retarder 13 will also 15 
work, as long as the two variable retarders oriented at differ- 
ent angles, however may involve more difficult calculations. 

In addition to the simplified measurement algorithm and 
auto-calibration, a key component in this polarization imag- 
ing apparatus is the variable phase retarder. The other means 0(| 
to change the phase retardation, rotating waveplates, Faraday 
rotators and Pockel’s cell have been used in laboratories 
worldwide. These devices, however, are bulky, high power 
consumption, slow or expensive. 

A preferred variable phase retarder is made from an elec- 
tro-optic (EO) material. The general requirement for an EO 
variable phase retarder is that, when a voltage is applied, a 
predetermined phase shift is produced for the light beam. 
Also the material is isotropic with no voltage applied, there is 
no residual birefringence to cause the initial phase shift. Pref- , f) 
erably, the material has a high electro-optic coefficient in 
order to reduce operating voltages to less than 500 volts. Also 
the mechanical characteristics allow formation of a bar or 
plate for use as the electrode. Of course, the material must be 
transparent at the wavelength of interest, e.g., between 450 , 
nm and 7 pm for most of the polarization imaging applica- 
tions. 

These requirements are satisfied by a class of ferroelectric 
complex oxides which 1) are optically isotropic; 2) have a 
Curie temperature less than about 490° C., so that electro- 4Q 
optic coefficients are high near room temperature; 3) have a 
diffusive phase transition, so that the temperature dependence 
of the electro-optic coefficients is lessened; and 4) which are 
not permanently poled by moderate electric fields since mate- 
rials with a low Curie temperature that become permanently 45 
poled are less stable. Example material systems include elec- 
tro-optic lanthanum modified PZT (PLZT), lead manganese 
niobate (PMN), a solid solution of lead zinc niobate and lead 
tantalate (PZN-PT), and a solid solution of lead manganese 
niobate and lead tantalate (PMN-PT). Besides PLZT, without 
being an exhaustive list the following materials and their solid 
solution with or without lanthanum modification, may be 
used: Pb(Zr,Ti)0 3 , Pb(Mg,Nb)0 3 , PbTi0 3 , Pb(Zn,Nb)0 3 and 
PbTa0 3 . These materials are available from Coming Inc. 
(Corning, N.Y.) and Boston Applied Technologies, Inc. .. 
(Woburn, Mass.), referenced as OptoCeramic materials. 
These materials are electrical controllable, solid-state, high 
speed, low insertion loss, and low fabrication cost. 

A preferred OptoC eramic material, PMN-PT, for use in the 
present invention may be described by the general formula 

Pb 1 -AMtMg^Nbi^Ti I , ! ! :/./h , l /2 3xI /2 17 /Si i.U/z/K 

wherein x is between about 0.6 and about 0.95, z is between 
about 0 and about 0.08, and t is between about 0.30 and about 
0.36. In especially preferred electro-optic ceramic materials 65 
of the present invention, z is between about 0.01 and about 
0.06. In preferred electro -optic materials of the present inven- 
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tion, t may be between about 0.32 and about 0.34, and x may 
be between about 0.65 and about 0.90. 

According to the present invention, the preferred OptoCe- 
ramic PMN-PT is lanthanum modified solid solution of lead 
manganese niobate and lead tantalate (La:PMN-PT, or 
PLMNT) with a nominal 3.5/75/25 La/PMN/PT composi- 
tion. For reference purposes nominal values for this PLMNT 
are n=2.48 and R=1 1 .5xl0 -16 nr/V 2 at 1 .06 pm and n=2.458 
and R=6.86xl0“ 16 nr/V 2 at 1 .55 pm, both at 23° C. PLMNT 
has a polycry stalline structure with crystal sizes ranging from 
about 5 to 20 microns. In the form of hot-pressed ceramics, it 
is optically isotropic and exhibits little birefringence with 
zero applied voltage. 

An electrode geometry of phase retarder which takes 
advantage of these materials is illustrated by the transverse 
field configuration illustrated in FIG. 3 . A plate 30 , for use as 
the variable phase retarder 12 and 13 of FIG. 1, is shown 
which has thin metallized electrodes 31 and 32 on the both 
sides, respectively, of a block section 33 of EO material. The 
electrodes 31 and 32 are shown connected to a driving voltage 
34 provide by controller 102 shown in FIG. 1 . The electric 
field (designated E) is 90° to the direction of light propaga- 
tion. Since the effect is based on the electric field in the 
material, it is desirable to arrange the electrodes as close 
together as possible to minimize the control voltage. The 
optical axis of the EO element is perpendicular to the electric 
field. 

As an example, an EO phase retarder is made from a 
PLMNT plate with metallized electrodes, which lias dimen- 
sions where the width (w) between the electrodes is about 5 
mm, a thickness (t) of about 4 mm, and a length (1) of about 5 
mm. Such a dimensioned PLMNT plate has an insertion loss 
of about 0.1 dB and half-wave voltage of 320V for a 700 nm 
wavelength laser at 23° C. 

It will be understood that the particular geometry described 
is an example and that other device geometries can be used. 
Different dimensions of w and t can be selected to provide the 
same it phase shift with different control voltage. The other 
preferred PLMNT plate has dimensions of, but not limited to 
a width (w) of about 0.5-5 mm, a thickness (t) of about 1 .30-5 
mm, and a length (1) of about 1 .0-5.0 mm. 

Another preferred OptoCeramic material is PLZT with a 
nominal 9.5/65/35 La/Pb/Zr composition. This composition 
is known to be transparent in a range from 450 nm to 7 pm. 

Yet another preferred OptoCeramic material, PZN-PT, for 
use in the present invention may be described by the general 
formula 

Pb , y ,B;t > L:i,[('Zn i Nb , _ t ) x Tl , _ x ] , - z /40 3 +x/ 2 - 3 xt /2- xz /&+ 

3xtz/8 

wherein x is between about 0.5 and about 0.9, y is between 
about 0.05 and about 0.5, z is between about 0 and about 0.05, 
and t is between about 0.30 and about 0.36. In especially 
preferred electro-optic ceramic materials of the present 
invention, x is between about 0.65 and about 0.85, y is 
between about 0. 1 and about 0.2, z is between about 0.02 and 
0.04, and t is between about 0.32 and about 0.34. 

Referring to FIG. 4 is the measured variable phase retarder 
response. The optical response 41 to the applied control volt- 
age 42 is less than 1 microsecond. This means that the speed 
of setting and resetting of the VPR is as fast as 1 .0 microsec- 
ond. According to the present invention, the imaging appara- 
tus only needs 2 and 3 settings for the first and the second 
variable phase retarders 12 and 13 to complete a measurement 
cycle, respectively. The polarization sampling time could be 
less than 3 microseconds. This equals the system can generate 
a set of 4 Stokes component images at a speed of up to 300K 
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frames per second (fps). With a fast enough imaging sensor, a 
real-time Stokes polarization imaging system can be realized. 

In the present invention, another preferred EO material for 
variable phase retarders 12 and 13 are made from liquid 
crystals for a polarization imaging system working between 5 
the optical wavelengths 350 nm to 2000 nm. Using the 
invented control algorithm illustrated in FIG. 2, and relation- 
ship between sequential intensity images and Stokes images 
listed in Table 1 , a fast polarization imaging system using LC 
variable retarders can also be configured. to 

A polarization imaging apparatus shown in FIG. 1 is basi- 
cally a passive imaging system which using ambient light 
source to illuminate the sample. Referring to FIG. 5 is a 
schematic view of a preferred embodiment of an active polar- 
ization imaging apparatus in accordance with the present 15 
invention. It consists of an illumination unit 50 and an imag- 
ing sensing unit 500. The said illumination unit and imaging 
sensing unit are aligned in a reflection configuration. Both 
illumination unit 50 and imaging sensing unit 100 can be set 
to an angle related to sample 56 surface normal 57. In addition 20 
to a light source 51, the illuminator may consist of parts such 
as a collimator 52, a polarizer or polarization modulator 53, 
and a fixed wavelength filter or tunable wavelength filter 54. 
Tunable optical band-pass filter was used to produce desired 
wavelength for the imaging, from visible to infrared (IR). The 25 
illuminator could produce a linear or arbitrary polarized light 
beam. 

A preferred light source 51 is a broadband light source, 
such as a tungsten-halogen lamp which provides a strong 
intensity over a broad spectrum from UV to NIR (300 nm to 30 
2 pm). Another preferred light source is a light emitting diode 
(LED) source or laser source. 

The scattering light 58 from sample 56 is detected by the 
imaging sensing unit 500. The optical assembly 11 in 100 
shown in FIG. 1 can be an optical objective lens 59 in FIG. 5 35 
with an infinite focal length, which makes the polarization 
imaging system a Stokes microscope. 

The first VPR 12 and second VPR 13 were controlled 
independently by a computer 102 through a controller unit 
101 which generates a sequential of voltages to control the 40 
phase retardations ofVPRs 12 and 13. A set of four intensity 
images, I 0 , I 1; I 2 and I 3 of the sample were captured by imag- 
ing sensor 15 when the phase retardations ofVPRs 12 and 13 
were set at (0,0), (it,0), (jt,Jt) and (jt/2,jt), respectively Then 
four Stokes components of a Stokes image, S 0 , Sj, S 2 and S 3 45 
can be calculated by computer 102 according to equation (1). 

Referring to FIG. 6 is a schematic view of a polarization 
imaging apparatus shown in FIG. 5 configured in a transmis- 
sion configuration. 

It will be apparent to those skilled in the art that various 50 
modifications and variations can be made to the present 
invention without departing from the spirit and scope of the 
invention. Thus it is intended that the present invention cover 
the modifications and variations of this invention provided 
they come within the scope of the appended claims and their 55 
equivalents. 

What is claimed is: 

1. A polarization imaging apparatus for measuring polar- 
ization images of a sample comprising: 

an optical assembly; 60 

a first variable phase retarder having a first controllable 
phase retardation; 

a second variable phase retarder optically coupled to the 
first variable phase retarder, the second variable phase 
retarder having a second controllable phase retardation; 65 
a linear polarizer optically coupled to the second variable 
phase retarder; 


an imaging sensor optically coupled to the polarizer for 
sensing the intensity images of the sample through the 
linear polarizer; 

a controller for varying the first controllable phase retarda- 
tion of the first variable phase retarder and for varying 
the second controllable phase retardation of the second 
variable phase retarders; and 

a computer for controlling the controller, recording the 
sequential intensity images, calibrating the alignment 
errors of the first and second variable phase retarders, 
and converting four sequential intensity images to four 
components of a Stokes image, 
wherein the computer is configured to calibrate the align- 
ment errors of first and second phase retarders using one 
known input polarization status and a set of sequential 
intensity images when setting the first and second vari- 
able phase retarders at different pairs of phase retarda- 
tions. 

2. The known input polarization status of claim 1 is a linear 
polarization. 

3. The known input polarization status of claim 1 is a 
circular polarization. 

4. The known input polarization status of claim 1 is an 
elliptical polarization. 

5. A set of sequential intensity images of claim 1 wherein 
the images are measured while the first and second variable 
phase retarders having 5 different sets of phase retardations. 

6 . A set of sequential intensity images of claim 1 wherein 
the images are measured while the first and second variable 
phase retarders having phase retardations of (jt/2,0), (jt,0), 
(O.jt), (jr.jt) and (ji/2,ji), respectively. 

7. A method to measure Stokes images of a sample, the 
method comprising the steps of: 

orienting a first variable phase retarder approximately at a 
22.5-degree angle with respect to the linear polarizer; 
orienting a second variable phase retarder approximately at 
a 45-degree angle with respect to the linear polarizer; 
illuminating the sample with light having a known polar- 
ization status; 

sensing a sequential of five intensity images with five sets 
of phase retardations while varying the first controllable 
phase retardation of the first variable phase retarder and 
varying the second controllable phase retardation of the 
second variable phase retarders; 
calculating the optical alignment errors of the first variable 
phase retarder and the second variable phase retarders; 
sensing a sequential of four intensity images with four sets 
of phase retardations while varying the first controllable 
phase retardation of the first variable phase retarder and 
varying the second controllable phase retardation of the 
second variable phase retarders; and 
converting a sequential of four intensity images to four 
components of a Stokes image. 

8. The method in claim 7 wherein the sample is illuminated 
by a linear polarized light. 

9. The method in claim 7 wherein the sample is illuminated 
by a circular polarized light. 

10. The method in claim 7 wherein the sample is illumi- 
nated by an elliptical polarized light. 

11. The method in claim 7, wherein alignment errors of first 
and second phase retarders are calibrated using the five sets of 
phase retardations of first and second variable phase retarders 
having the phase retardations of (jc/ 2,0), (jt,0), (0 ,jt), (it, it) and 
(jt/2,jr), respectively. 
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12. The method in claim 7 further comprising the steps of: 
illuminating the sample with a light source; 
calibrating the ali gnm ent error and phase retardation error 
of the phase retarders; 

adjusting the light source polarization state; 5 

changing the light source wavelength spectrum; 
acquiring the intensity images; and 
converting to Stokes images. 



